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For the hydrogen atom in combined magnetic and electric 
fields we investigate the dependence of the quantum spectra, 
classical dynamics, and statistical distributions of energy lev- 
els on the mutual orientation of the two external fields. Reso- 
nance energies and oscillator strengths are obtained by exact 
diagonalization of the Hamiltonian in a complete basis set, 
even far above the ionization threshold. At high excitation 
energies around the Stark saddle point the eigenenergies ex- 
hibit strong level repulsions when the angle between the fields 
is varied. The large avoided crossings occur between states 
with the same approximately conserved principal quantum 
number, n, and this intramanifold mixing of states cannot be 
explained, not even qualitatively, by conventional perturba- 
tion theory. However, it is well reproduced by an extended 
perturbation theory which takes into account all couplings be- 
tween the angular momentum and Runge-Lenz vector. The 
large avoided crossings are interpreted as a quantum mani- 
festation of classical intramanifold chaos. This interpretation 
is supported by both classical Poincare surfaces of section, 
which reveal a mixed regular-chaotic intramanifold dynamics, 
and the statistical analysis of nearest-neighbor-spacing distri- 
butions. 

PACS numbers: 32.60. -fi, 03.65. Sq, 05.45. -Kb, 32.70.Cs 



I. INTRODUCTION 

The effects of external electric and magnetic fields 
on atomic spectra is a fundamental question of atomic 
physics. However, up to now most investigations have 
concentrated on atoms in only one of the two fields, or on 
selected mutual orientations such as parallel or perpen- 
dicular fields. For example, Rydberg atoms in magnetic 
fields are nonintegrable systems and have been shown to 
be ideally suited for studies of the quantum manifesta- 
tions of classical chaos [0-^ . Because of the cylindrical 
symmetry around the magnetic field axis, one component 
of the angular momentum is conserved, and the problem 
is nonseparable in two degrees of freedom. The cylin- 
drical symmetry is broken in combined nonparallel mag- 
netic and electric fields, and the system becomes non- 
separable in three degrees of freedom. The special case 
of perpendicular fields has been investigated both exper- 
imentally and theoretically by quantum |^-|^ and 
classical [[lO[-p3[ methods. However, the case of arbitrary 
field orientation is the most general situation for atoms in 
uniform external fields, and therefore deserves appropri- 
ate attention. For the general field arrangement quantum 



calculations have been performed so far for weak external 
fields in first and second order perturbation theory 
and in the regime of very strong magnetic and electric 
fields 

In this Paper we investigate the hydrogen atom in ex- 
ternal fields with arbitrary mutual orientations. In Sec. II 
we calculate the (numerically) exact quantum mechanical 
eigenenergies and oscillator strengths for a wide energy 
and field range, even far above the ionization threshold 
(Stark saddle point), by diagonalization of the Hamilto- 
nian matrix in a complete set of eigenfunctions. On the 
basis of these exact quantum calculations we investigate 
the dependence of resonance energies on the angle be- 
tween the external fields. At energies around the ioniza- 
tion threshold we observe large avoided crossings, which 
cannot be explained even qualitatively by conventional 
perturbation theory. These avoided crossings mainly oc- 
cur between eigenstates of the same approximately con- 
served principal quantum number n. The phenomenon 
is interpreted as a quantum manifestation of intramani- 
fold chaos, which has recently been discovered in classi- 
cal investigations of the hydrogen atom in perpendicular 
crossed magnetic and electric fields pT[ |. 

To verify the interpretation of the avoided crossings 
as intramanifold chaos we extend, in Sec. Ill, the con- 
ventional perturbation theory by taking into account 
all couplings between the angular momentum and the 
Runge-Lenz vector. This extended quantum mechanical 
perturbation theory can well reproduce, at least qualita- 
tively, the large avoided crossings between levels of the 
same principal quantum number n. 

Based on the Hamiltonian of the extended perturba- 
tion theory we analyze, in Sec. IV, the classical motion 
of the angular momentum and the Runge-Lenz vector by 
means of Poincare surfaces of section (PSOS). The PSOS 
reveal large chaotic fractions of the classical phase space 
for those energy-field regions and mutual field orienta- 
tions where the quantum spectra exhibit large avoided 
crossings. The classical calculations therefore confirm the 
interpretation of the intra n-manifold level repulsion as 
a quantum manifestation of intramanifold chaos. 

In Sec. V we investigate the nearest-neighbor-spacing 
distributions of states with the same principal quantum 
number n. The distributions show a Poissonian behav- 
ior for parallel fields, where the intramanifold dynam- 
ics is regular. For arbitrarily oriented fields the distri- 
bution turns to a Brody-type distribution, indicating a 
mixed regular-chaotic intramanifold dynamics, in agree- 
ment with the classical Poincare surface of section anal- 
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ysis. 



II. EXACT QUANTUM CALCULATIONS 

The electronic motion of the hydrogen atom in uniform 
magnetic and electric fields with arbitrary mutual orien- 
tations is nonseparable in three degrees of freedom. The 
Hamiltonian [in atomic units, / — F/(5.14 x 10^ V/cm), 
7 = B/{2.ib X 10^ T)] reads 

= ^P' - ^ + \lL, + \-i\x^ + y2) + /^a: + , 

(1) 

where the magnetic field is oriented in the z-direction and 
the electric field in the {x, z)-plane with /|| and /j^ the 
components of the electric field parallel and perpendic- 
ular to the magnetic field axis, respectively. We neglect 
relativistic and spin effects as well as effects due to the 
finite nuclear mass, which yield only very small contribu- 
tions in the energy-field regions we examine. 

For the exact quantum calculations of the hydrogen 
atom in external fields with arbitrary mutual orientations 
we extend the method described in Ref. Q for the special 
case of the hydrogen atom in perpendicular magnetic and 
electric fields (/|| =0). In perpendicular fields the parity 
with respect to the {z = 0)-plane, tt^, is an exact quan- 
tum number and allows diagonalizing the Hamiltonian 
in each subspace tt^ = separately. This symmetry 
is broken by the parallel component of the electric field 
resulting in an increase of the basis size for the numerical 
diagonalizations by about a factor of 2. 

For the numerical calculations we transform the Hamil- 
tonian by introducing dilated semiparabolic coordinates 
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where 6 is a free length scale parameter. The Schrodinger 
equation in dilated semiparabolic coordinates then reads 
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The Schrodinger equation (||) has the form of two cou- 
pled two-dimensional harmonic oscillators, and is there- 
fore conveniently represented, in matrix form, in terms of 
the complete set of basis functions given by the product of 
the eigenstates \Npm) of the two-dimensional harmonic 
oscillator, viz. 



\N^,N^m) = \N^m) x \N^m) 



(6) 



One is led to a generalized eigenvalue problem with sparse 
symmetric matrices. The elements of these matrices can 
be calculated using the familiar operator relations for the 
two-dimensional harmonic oscillator. More details of the 
method are described in p^ ]. 

To account for continuum states we adopt the complex- 
rotation method [ pT|Jl^ , which is based on the replace- 
ment of the coordinate vector r with re*^ in the Hamil- 
tonian and the wave function, and has proved very effi- 
cient in the calculation of resonances above the ionization 
threshold for the hydrogen atom in magnetic fields |l9| ] 
and in crossed magnetic and electric fields ||]. By this 
transformation, hidden resonances of the Hamiltonian in 
the continuum, associated with complex eigenvalues, are 
exposed, while the resonance wave functions can still be 
described by the £2 integrable basis functions (^), but 
with complex arguments. In our approach, the complex 
rotation by the angle 9 corresponds to the complex di- 
latation 
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in (H) and in Schrodinger's equation (||). 

Numerically a generalized eigenvalue problem with 
complex symmetric non-Hermitian matrices has to be 
solved, which was achieved by extending the spectral 
transformation Lanczos method [po[ to complex matri- 
ces. From the imaginary part of the complex energies, 
E, we obtain the corresponding widths, F, and lifetimes, 
T, of the resonances by 



T 



= -2 Im E 



(8) 



As an example we calculated spectra of the hydro- 
gen atom in combined magnetic and electric fields with 
B = 100 T, F = 50 kV/cm. These field strengths are 
by about an order of magnitude larger than typical labo- 
ratory fields. However, they allow studying the interest- 
ing physical effects at energy regions with a reasonably 
low density of states and reduce the computer resources 
(CPU time and storage requirements) necessary for the 
numerical calculation of the spectra. Bound states and 
resonances were obtained by numerical diagonalization 
of Eq. I in a basis set of dimension < 6201. The com- 
plex energies are presented in Fig. 1 for an angle /3 = 
arctan(/||//x) = 7r/4 between the fields. The energy of 

the Stark saddle point, Esp = —2^y\f \ au = —1368 cm^^ 
is marked by an arrow in Fig. 1 . Below the threshold all 
eigenvalues are located on the real axis (see Fig. 1). 
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FIG. 1. Eigenenergies and resonances in the complex en- 
ergy plane for the hydrogen atom in magnetic and electric 
fields {B = 100 T, F = 50 kV/cm) of mutual orientation 
/3 = 45°. The energy of the Stark saddle point is marked by 
an arrow. 

Above threshold, long-lived states close to the real en- 
ergy axis still exist, but hidden resonances of the Hamil- 
tonian, associated with complex eigenvalues, are exposed 
by the complex rotation method. 

With the energy eigenvalues and eigenvectors of (||) at 
hand it is also possible to calculate the cross section for 
dipole transitions from an initial state ^'o with energy 
Eq. The cross section can be written in the form pit] 
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FIG. 2. (a) Photoabsorption spectrum for the hydro- 
gen atom in magnetic and electric fields {B = 100 T, 
F = 50 kV/cm) of mutual orientation (3 = 45°. Transitions 
\2p0) IV'/) with light polarized parallel to the magnetic 
field axis, (b) Resonances in the complex energy plane. 
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where j{9) are final states at complex energies Ej{6), D 
denotes the dipole operator for some given polarization, 
and a « 1/137 is the fine-structure constant. As an 
example a photoabsorption spectrum for transitions from 
the initial state \2p0) with light polarized parallel to the 
magnetic field axis is presented in Fig. 2a. The chosen 
energy region is the high energy part of Fig. 1 and is well 
above the Stark saddle point, Esp — —1368 cm"-'^. The 
spectrum clearly exhibits sharp peaks of long-lived states 
superimposed on broad line shapes of rapidly decaying 
resonances and a continuous background. The complex 
energy eigenvalues are shown in Fig. 2b for comparison. 

For atoms in external fields with arbitrary mutual ori- 
entations the angle /3 between the magnetic and electric 
field axis is an additional free parameter which allows 
one to continuously vary the geometry between the ex- 
treme situations of parallel and perpendicular fields. It 
is therefore interesting to fix the absolute values of the 
electric and magnetic field strength and to follow the en- 
ergy levels as a function of the angle (3. A part of such an 
{E, ^)-diagram is shown in Fig. 3 at B = 100 T, = 50 
kV/cm. The dotted line marks the energy of the Stark 
saddle point. The (i?, /3)-diagram exhibits quite com- 
plicated patterns which cannot be explained, not even 
qualitatively, by conventional perturbation theory B . 
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FIG. 3. Energy-angle diagram for the hydrogen atom in 
magnetic and electric fields at B = 100 T, _F = 50 kV/cm. /3 
is the angle between the two field axis. 

In particular, many levels undergo large avoided cross- 
ings, which are found to be most pronounced around an- 
gles 40° < (3 < 70°. The avoided crossings indicate chaos 
in the underlying classical dynamics. It is the purpose of 
the following Sections to investigate and analyze the pat- 
terns observed in Fig. 3 in more detail. This is performed 
by means of an extended quantum and classical pertur- 
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bation theory, which will allow the interpretation of the 
avoided crossings as effects of intramanifold chaos. 



III. EXTENDED QUANTUM MECHANICAL 
PERTURBATION THEORY 

The perturbation theory of the hydrogen atom in com- 
bined electric and magnetic fields goes back to the early 
days of quantum mechanics. The first order perturba- 
tion theory was developed by Epstein Born [ p3[ , 
and Pauli 0. The second order equations describing 
the eigenenergies up to the quadratic terms in the field 
strengths were derived by Solov'ev The conventional 
perturbation theory is based classically on the construc- 
tion of adiabatic constants of motion which are related 
to approximate quantum numbers in quantum mechan- 
ics. The results of second order perturbation theory agree 
well with our exact quantum calculations for weak exter- 
nal fields and at low excitation energies. However, at 
strong perturbations, e.g., at excitation energies around 
the Stark saddle point, the conventional perturbation 
theory cannot explain the breakdown of approximate 
quantum numbers and the onset of chaos in the classical 
dynamics. Therefore the large avoided crossings which 
have been observed, e.g., in Fig. 3 are not reproduced by 
the conventional perturbation theory. In the following we 
will briefiy review the conventional perturbation theory 
and then generalize the method by taking into account 
all couplings between states within the same n-manifold 
up to second order in the external fields. The concept 
does not assume the existence of a complete set of ap- 
proximate quantum numbers and therefore can describe 
and explain the phenomenon of intramanifold chaos. 

The perturbation series for the hydrogen atom in ex- 
ternal magnetic and electric fields up to second order in 
the field strengths reads 



and the Runge-Lenz vector 



Hn^-T^+Vi+V2 + W 



with 



Vi = -jL, + f • r 
the paramagnetic and linear Stark terms, 

V2 = l^f^P^ 
the diamagnetic term, and 

1^ = (f •r)G„(f -r) 



(10) 



(11) 



(12) 



(13) 



the quadratic Stark effect, with G„ the Green's function 
of the field free hydrogen atom, Hq. The operators Vi, 
V2, and W can now be replaced within each n-manifold 
by operator identities in terms of the angular momentum 



1 



1 r 
-(pxL-Lxp)-- 

2 r 



or in terms of the linear combinations of these vectors 



(15) 



Ii,2 = i(L±A) 



(16) 



The vectors Ii and I2 commute with each other and their 
components obey the commutation relations for angular 
momentum operators. 



[h,j2k] = 0, 

[/ij, life] = iejkilii , 
[l2j,hk] — if-jkihi ■ 



(17) 



Both vectors have the same norm, which depends on the 
principal quantum number. 



(18) 



The operator identity for Vi in a given n-manifold then 
reads M 



with 



Vi uji ■ li + l32 ■ I2 



wi,2 = 2 (7 T 3n/ 



(19) 



(20) 



as is illustrated in Fig. 4. Because Ii and I2 indepen- 
dently fulfill the commutation algebra of two angular 
momenta, the operator Vi in (|l9|) can be quantized im- 
mediately to obtain the energy correction in first order 
perturbation theory 



E^^\ „ = com' 

nn'n" ^ 



- ijJ2n 



(21) 



with uii^2 = |wi.2| and quantum numbers n', n" — — (n — 
l)/2, -(n - 3)/2, . . . , +(n - l)/2. Classically, Ii and I2 
turn around the vectors wi and ijJ2 , whose absolute values 
are the frequencies of the secular motion (see Fig. 4). 

The operator identities for V2 and W are the starting 
point for the second order_perturbation theory and have 
been derived by Solov'ev 



V2 



16 
16 



5 + + _ 3Al) , (22) 
(Sn^ 3H- 24L2 - 21L^ + 9A}) , (23) 



with Lf and Af the projections of the angular momen- 
tum and Runge-Lenz vector on the electric field axis. 

Except for the special case of perpendicular fields 
where we have wi = W2 = 5\/7^ + Sn^/^, the degen- 
eracy of the states belonging to each n-manifold is 



L = r X p 



(14) already completely destroyed in first order perturbation 
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FIG. 4. Schematic view of the vectors lui, 102 and Ii, I2 
which are used in perturbation theory. 

theory. In conventional perturbation theory, the second 

(2) 

order energy correction, E^^,^,, is given as the expecta- 
tion value of V2 

^(2) 

'nn' 

= j^f!''"^ + 19-12 [n''^ + n'n" cos{ai + aa) 

Tl^ 2 

+ n""^)] H [1n^ + 5 + Ann" sin ai sin Qf2 

48 

+ {n^ - l)(cos^ ai + cos^ 02) - 12 (n'^ cos^ ai 

— n'n" cos ai cos 02 + n"'^ cos^ 0:2)] , (24) 

where ai and 0:2 are the angles between the magnetic 
field axis and the vectors uji and 0J2, respectively (see 
Fig. 4). The energy eigenvalues are finally obtained as 

Enn'n" - -l/27i2 + + E^^' n" ■ Notc that Eq. 

is not valid for perpendicular fields because in this 
case the degeneracy of the n-manifolds is not completely 
destroyed in first order perturbation theory How- 
ever, the qualitative results of the following discussion are 
not changed when the conventional perturbation theory 
for nearly perpendicular fields is applied instead of Eq. 

HQ- 

The results of the conventional second order perturba- 
tion theory for non-perpendicular fields are compared to 
the exact quantum calculations in Fig. 5. The dashed 
lines in Figs. 5a and 5b are the levels of manifolds n = 9 
and n — 10, respectively. To simplify the identification 
the exact levels with the corresponding principal quan- 
tum numbers are highlighted by solid lines while other 
levels are drawn by dotted lines. Although the identi- 
fication is possible, the comparison reveals remarkable 
discrepancies between the results of conventional pertur- 
bation theory and the exact quantum results in this en- 
ergy regime. In the conventional perturbation theory all 
level crossings are exact and cannot reproduce the large 
avoided crossings discussed above. Another interesting 
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FIG. 5. Exact energy eigenvalues (solid lines) compared 
with energy values calculated by conventional second order 
perturbation theory (dashed lines), (a) Principal quantum 
number n = 9; (b) n — 10. 

result of the comparisons in Fig. 5 is that the large 
avoided crossings do not occur between levels with dif- 
ferent principal quantum numbers. Therefore they can- 
not be interpreted as an n-mixing effect, which has been 
identified as the origin of level repulsion in the hydro- 
gen atom in magnetic or parallel magnetic and electric 
fields. The level repulsions in Fig. 5 solely occur between 
levels of the same n-manifold and therefore must be inter- 
preted as a mixing between states \nn'n") with different 
quantum numbers n' and n" but with the same princi- 
pal quantum number n. The breakdown of the quantum 
numbers n' and n" is classically related to intramanijold 
chaos, which will be discussed in Sec. IV. 

We now demonstrate that the avoided crossings are an 
effect of intramanifold mixing of states. The conventional 
perturbation theory is based on the assumption that the 
external fields are weak, i.e., the linear terms in the ex- 
ternal fields, Vi, and the quadratic terms, V2 -I- W , can 
be treated separately as first and second order pertur- 
bations. We now extend the conventional perturbation 
theory by diagonalizing the complete perturbation oper- 
ator Vi + V2 + 14^ in Eq. |l^ within a given n-manifold. As 
a basis set we choose the parabolic states, i.e., simulta- 
neous eigenstates of the Coulomb Hamiltonian, Hq, and 
the z components of the angular momentum and Runge- 
Lenz vectors, and A^. By substituting the angular 
momentum and Runge-Lenz vector in Eqs. |l^, |2^, and 
E3 with L = Ii + I2 and A = Ii — I2 and applying 
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the commutator algebra for the operators Ii and I2 it is 
a straightforward task to compute the required matrix 
elements oiVi + V2 + W. The energy eigenvalues are ob- 
tained by numerical diagonalization of the perturbation 
matrix of dimension x v? . 

The results of the extended perturbation theory are 
presented in Fig. 6 for parts of the manifolds n — % and 
n = 10. The exact quantum levels (highlighted solid and 
dotted lines) are the same as in Fig. 5. 
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FIG. 6. Same as Fig. 5 but for the comparison between the 
exact eigenenergies (solid lines) and the results of extended 
perturbation theory (dashed lines), (a) n = 9; (b) n = 10. 

The dashed lines are the eigenvalues from the diagonal- 
ization of the complete perturbation matrix ([l^) . They 
do not totally agree with the exact levels, but this can- 
not be expected at the high excitation energies close to 
the Stark saddle point. However, Fig. 6 qualitatively ex- 
hibits the same level dynamics for both the solid and 
dashed lines for all mutual field orientations, i.e., from 
parallel fields (/3 = 0) to perpendicular fields {(3 — 7r/2). 
In particular, the avoided crossings between levels of the 
same n-manifold are well reproduced within the extended 
perturbation theory, and this strongly supports the inter- 
pretation of the observed level repulsion as an effect of 
intramanifold mixing. The n-mixing of states will be- 
come important only at higher energies far above the 
ionization threshold. 



IV. CLASSICAL INTRAMANIFOLD DYNAMICS 

The hydrogen atom in magnetic and electric fields with 
arbitrary mutual orientations is nonintegrablc in three 



degrees of freedom. However, for weak external fields 
and low excitation energies three approximate quantum 
numbers n, n' , and n" exist. In classical mechanics these 
quantum numbers are related to three constants of the 
motion, viz. the length of the vectors Ii and I2, 

|Ii| = |l2| = ^n (25) 

(note the slight difference to the quantum expression 
(|l8|)) and their projections on the axis ui and L02 defined 
in Eq. |2^. Ii and I2 turn around the vectors di and ijJ2 
with constant frequencies (see Fig. 4), i.e., the classical 
dynamics is completely regular. In the extended pertur- 
bation theory discussed in the previous Section, states 
\nn'n") with the same principle quantum number n but 
with different quantum numbers n' and n" are mixed, 
i.e., two approximate quantum numbers are destroyed. 
The breakdown of two constants of motion can lead to a 
chaotic secular motion of the vectors Ii and I2 even when 
the principal quantum number, n, is still conserved. This 
phenomenon is called intramanifold chaos and has been 
discovered in the hydrogen atom in perpendicular crossed 
magnetic and electric fields . Here we study the clas- 
sical intramanifold dynamics of the hydrogen atom in ex- 
ternal fields with arbitrary mutual orientations, and are 
especially interested in the dependence of the classical 
dynamics on the angle between the external fields. The 
quantum calculations in Sec. Ill revealed large avoided 
crossings at angles 40° < (i < 70°. If the level repul- 
sion is the quantum counterpart of classical intramani- 
fold chaos, the classical dynamics should be particularly 
irregular at those field orientations. 

To investigate the intramanifold dynamics we solve the 
equations of motion for the secular motion of the two vec- 
tors Ii and I2 defined in (^|). By expressing the angular 
momentum and Rungc-Lcnz vector in the Hamiltonian 
( |lo| ) in terms of Ii and I2 , scaled with the principal quan- 
tum number, viz. Ii,2 — > Ii,2/"-, we obtain the classical 
Hamiltonian valid within a given n-manifold, 

= -1 + n^jihz + hz) 

- 3n^/[sin/3(/ij; - hx) + cos/3(/i2 - l2z)] 

+ ^("'7)' [3 - 4.(11 + II - hzhz) - Shxhx - &hyl2y] 
O 

- i (nV)'{17 - 12 sin^ Pill + II) - 12 cos' f3{ll + /|J 

+ (48 - 60 sin^ I3)hxhx + -iShyhy + (48 - 60 cos^ p)hzl2z 

- 12sin/3cos/3[5(/i^/2z + hzhx) + 2{h Jiz + hxhx)]} ■ 

(26) 

The classical equations of motions are now obtained from 
the Hamiltonian (|2^) by the Poisson brackets 

11 = , 

12 = [l2,H], (27) 
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which can be solved with the help of the elementary Pois- 
son brackets for the components of Ii and I2, 



[hj,hk\ = 0, 
[hjihk] — ^jkihi 



(28) 



The classical intramanifold dynamics depends on the en- 
ergy and field strengths scaled with the principal quan- 
tum number, ri^E, n^j, and n*/, and on the angle, be- 
tween the fields. The classical trajectories are computed 
in the six-dimensional space {Ii(^), l2(t)}. However, the 
norm of Ii and I2 is conserved (pq), and therefore the 
effective phase space is four-dimensional. A convenient 
parameterization of the four-dimensional phase space are 
the projections of Ii and I2 on the axis uji and UJ2 (see 
Fig. 4), called Jiz and J2Z in what follows, and the polar 
angles 0i and (f>2 of Ii and I2 with respect to the uji^2 
axis. Note that {Jiz,4'i} and {J2z,(t>2} are action-angle 
variables and Jiz and J22 are constants of motion in the 
limit of weak non-perpendicular external fields. 

A common way of visualizing the classical dynamics 
of a system is the method of Poincare surfaces of section 
(PSOS). In a system with n degrees of freedom the PSOS 
is in general a (2n — 2)-dimensional subspace of the (2n)- 
dimensional phase space. Elliptic fixpoints surrounded 
by torus structures indicate regular classical dynamics 
while chaotic motion is indicated by stochastic layers in 
the PSOS. For the analysis of the intramanifold dynam- 
ics we use the action-angle variable of the weak field limit 
and define the intersections with the PSOS by (j>2 — 0. 
The two-dimensional PSOS is now defined by the doublet 
{•/iz,0i}7 while J22 is given implicitly via the conserva- 
tion of energy. 
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As an example we calculated PSOS for very weak ex- 
ternal fields, — 0.053, n'^f — 0.0061, which is related 
to the manifold n = 5 at the field strengths i? = 100 
T, F = 50 kV/cm chosen for the quantum calculations 
in Sees. II and III. Figs. 7a and 7b present the results 
for mutual field orientations (3 = (parallel fields) and 
(3 = 7r/4, respectively. Qualitatively both PSOS look 
the same, i.e., they exhibit nearly parallel lines with 
Jiz ~ const. This is a verification of the fact that Jiz is 
a constant of motion in the limit of weak external non- 
parallel fields. For the special case of perpendicular fields 
(/? = ""/S) the PSOS have a qualitatively different ap- 
pearance (Fig. 7c) . This is a consequence of the fact that 
in perpendicular fields the degeneracy is not completely 
destroyed in first order perturbation theory (lui = lu2 in 
Eq. |2^), and therefore Jiz is no constant of motion. How- 
ever, Fig. 7 clearly reveals that the classical intramanifold 
dynamics is completely regular in the weak field limit for 
all mutual field orientations. 

The situation changes when we look at the PSOS for 
stronger interaction of the hydrogen atom with the ex- 
ternal fields. To search for classical intramanifold chaos 
we analyzed the dynamics at parameters n^j = 0.74, 
n^/ = 0.2, Ti^E = —0.5, which is related to the manifold 
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FIG. 7. Poincare surfaces of section for the classi- 
cal intramanifold dynamics with parameters = 0.053, 
n"/ = 0.0061, ri^E = -0.5. Mutual field orientation (a) 

(3 = 0°, (b) /3 = 45°, (c) (3 = 90°. 



FIG. 8. Same as Fig. 7 but with parameters 71^7 = 0.74, 
n^f = 0.2, n^E = —0.5 and mutual field orientation (a) 

(3 = 20°, (b) (3 = 58°, (c) (3 = 60°, (d) P = 62°, (e) (3 = 90°. 
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n = 12 at the field strengths 5 = 100 T, F = 50 kV/cm. 
The Poincare surfaces of section are presented in Fig. 8 
for various mutual orientations of the external fields. In 
the region < (3 < 40° the torus structures Jiz « const 
observed at weak fields (Figs. 7a and 7b) are strongly 
distorted (see Fig. 8a at /? = 20°). However, the tori are 
not destroyed, i.e., the classical intramanifold dynamics is 
mostly regular apart from tiny stochastic regions around 
hyperbolic fix points. Note that no intramanifold level 
repulsion was observed for those field orientations in the 
quantum calculations of Sees. II and III. 

The onset of strong classical intramanifold chaos is dis- 
covered at mutual field orientations 40° < (3 < 70°. Ex- 
amples at /3 = 58°, (3 — 60°, and /3 — 62°, are given 
in Figs. 8b, 8c, and 8d, respectively. They clearly ex- 
hibit mixed regular-chaotic classical dynamics with large 
or even dominating stochastic, i.e., chaotic phase space 
regions. Fig. 8 also shows a sensitive dependence of 
the classical dynamics on the angle between the external 
fields. Even small deviations A/3 = 2° of the mutual field 
orientation result in significant changes of the PSOS. For 
example at /3 = 58° (Fig. 8b) there exists a regular region 
at the top of the PSOS, which vanishes at (3 = 60° in Fig. 
8c. At that angle the PSOS exhibits two large irregular 
areas which are dynamically separated. The dynamical 
barrier is resolved in Fig. 8d where only a few stabil- 
ity islands are embedded in one large irregular region. 
The field arrangement where strong classical intraman- 
ifold chaos is observed (40° < (3 < 70°) coincides with 
the region where large avoided crossings between quan- 
tum levels with the same principal quantum number, n, 
were found in Fig. 6. This strongly supports the interpre- 
tation of the level repulsion as a quantum manifestation 
of intramanifold chaos. 

If the angle between the fields is further increased 
{[3 > 70°) the intramanifold dynamics becomes more 
and more regular again. The PSOS for the special case 
of perpendicular fields (/3 — 90°) is presented in Fig. 8e. 
Chaotic regions exist in agreement with previous classical 
investigations of the crossed field atom ||l^, however, the 
PSOS is clearly dominated by regular torus structures. 



V. INTRAMANIFOLD LEVEL STATISTICS 

It is well established that the nearest-neighbor-spacing 
distribution (NNS) of integrable quantum systems is 
given, after unfolding the spectra to unit mean level spac- 
ing (s) = 1, by a Poisson distribution 



Pi 



Poisson 



(s) = exp [-s] 



(29) 



while quantum systems with a fully chaotic (ergodic) 
underlying classical dynamics are characterized by the 
Wigner distribution 



obtained from random matrix theory |25 2^. In sys- 
tems with a mixed regular-chaotic classical dynamics 
the nearest-neighbor-spacing distribution can be phe- 
nomenologically described by the Brody distribution |2^] 



PBrody(s; q) = {q + l)as« exp (31) 



where 









[ V9H 





9+1 



(32) 



Avignor(s) = Y^^P 



(30) 



and q is a parameter which interpolates between the Pois- 
son distribution (g = 0) and the Wigner distribution 
(q = 1) and is roughly related to the percentage of chaotic 
phase space volume of the underlying classical system. 

Here we do not intend to compare nearest- 
neighbor-spacing distributions with the complete three- 
dimensional classical dynamics of the hydrogen atom in 
external magnetic and electric fields with arbitrary mu- 
tual orientations. Instead, we want to search for finger- 
prints of intramanifold quantum chaos in the level statis- 
tics. We therefore do not investigate the NNS distribu- 
tion of all levels of a given spectrum but analyze sepa- 
rately the spacings of neighboring states with the same 
principal quantum number, n. For direct comparisons 
with the PSOS of the classical intramanifold dynamics 
(see Sec. IV) the spacing distributions should be ob- 
tained at constant parameters n^/, and n'^E for 
the scaled field strengths and energy. In order to get a 
reliable statistics for the problem we have to calculate a 
sufficiently large set of eigenvalues. To fulfill these condi- 
tions we calculated eigenenergies at constant scaled field 
strengths and n*/ by means of the extended sec- 
ond order perturbation theory derived in Sec. HI. The 
Hamiltonian ( p^ was diagonalized numerically for all n- 
manifolds from n = 50 to n = 60. The resulting spectra 
were unfolded separately to mean level spacing (s) = 1 
psf . For the NNS distributions we included only eigen- 
statcs whose scaled energies deviate by at most O.lri'^7 
from the given value of n^E. This ensures that the scaled 
energy is approximately constant and the NNS distribu- 
tions can be compared to the classical dynamics, viz. the 
PSOS, at a given constant energy. 

As an example. Fig. 9 presents the NNS distribution at 
field parameters n^-f = 0.31, n'^f = 0.064, f3 = 45°, and 
for approximately constant scaled energy n^E « —0.5. 
These parameters represent the situation of the manifold 
n = 9 at B = 100 T and F = 50 kV/cm in Fig. 6a. How- 
ever, to increase the density of states for a more reliable 
statistics we analyzed the manifolds n = 50 to n = 60 
at correspondingly rescaled energies and field strengths. 
The histogram in Fig. 9 is close to a Poisson distribution 
( p9| ) indicating a regular intramanifold dynamics. An 
even better fit is obtained by a Brody distribution ( |3l| ) 
with q = 0.08 (dashed line in Fig. 9). The deviations 
from the Poisson distribution can be interpreted as the 
onset of intramanifold quantum chaos. The decrease of 
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FIG. 9. Intramanifold ncarest-neighbor-spacing distribu- 
tion for the parameters n^'y = 0.31, n* f = 0.064, ri^E ~ —0.5, 
(5 = 45°. The dashed hne is a Brody distribution with 
q = 0.08. 

the probability P{s) at small spacings s is related to the 
observation of level repulsions in Fig. 6a. For a more 
pronounced verification of intramanifold quantum chaos 
in NNS distributions we investigated quantum spectra 
at parameters n^7 = 0.74, n^f = 0.2, n^E « -0.5, i.e., 
the same parameters as for the analysis of the classical 
intramanifold dynamics in Fig. 8 of Sec. IV. The NNS dis- 
tributions for the same set of five different mutual field 
orientations between (3 = 20° and (3 = 90° are presented 
in Fig. 10. The histograms indicate a predominately reg- 
ular intramanifold dynamics at /? = 20° (Fig. 10a) and 
/3 = 90° (Fig. lOe) with Brody parameters g = 0.1 and 
q = 0.05, respectively, and this is in excellent agreement 
with the regular torus structures in the PSOS of Figs. 8a 
and Be. The NNS distribution changes dramatically at 
mutual field orientations around /3 = 60° (see Figs. 10b 
- lOd). The histograms turn into Brody distributions 
with much higher Brody parameters, up to g = 0.45, 
indicating a large irregular part of the underlying clas- 
sical intramanifold dynamics. Again the behavior of the 
NNS distributions is in very good agreement with the 
PSOS of Figs. 8b - 8d, which exhibit large stochastic 
regions at field orientations around (3 = 60°. Another 
interesting aspect which can be observed in Fig. 10 is the 
sensitive dependence of the NNS distributions on small 
variations of the mutual orientation between the external 
fields around /3 w 60°. The Brody parameter increases 
from g = 0.3 at /? = 58° (Fig. 10b) to g = 0.45 at /? = 60° 
(Fig. 10c) and then decreases to q = 0.2 ai j3 — 62° (Fig. 
lOd). The rapid changes of the Brody parameter reflect 
the sensitive dependence of the classical intramanifold 
dynamics on the field orientation discovered in the PSOS 
(see Figs. 8b - 8d). 



VI. CONCLUSION 

We have investigated the hydrogen atom in electric 
and magnetic fields with arbitrary mutual orientations 
by means of exact quantum calculations, an extended 
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FIG. 10. Intramanifold ncarest-neighbor-spacing distribu- 
tions for the parameters n^7 = 0.74, n* f = 0.2, n^E « —0.5 
and mutual field orientations (a) /3 = 20°, (b) /3 = 58°, (c) 
13 = m°, (d) 13 = &2°, (e) /3 = 90°. The dashed lines are 
fitted Brody distributions with (a) q = 0.1, (b) q = 0.2, (c) 
q = 0.45, (d) q = 0.3, (e) q = 0.05. 

second order perturbation theory, classical Poincare sur- 
face of section analysis of the intramanifold dynam- 
ics, and statistical analysis of the intramanifold nearest- 
neighbor-spacing distributions. For mutual orientations 
aroimd (3 w 60° and excitation energies aroimd the Stark 
saddle point the exact quantum calculation reveal large 
avoided crossings between states with the same approx- 
imate principal quantum mimber, n. Because the level 
repulsion occurs preferably within a given n-manifold it 
is interpreted as a quantum manifestation of intramani- 
fold chaos. This interpretation is strongly supported by 
the analysis of both the classical intramanifold dynamics 
and the intramanifold NNS distributions. The classical 
analysis exhibits large stochastic regions, i.e., chaotic mo- 
tion in the PSOS at those parameters where strong level 
repulsion has been observed in the quantum spectra, and 
the intramanifold NNS distributions show a transition 
from a Poisson distribution to a Brody distribution at 
those energy-field parameters. 

In parallel magnetic and electric fields {(3 = 0) the 
intramanifold dynamics is regular because of the cylin- 
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drical symmetry of the system and the conservation of 
the z-component of the angular momentum. It might be 
expected that intramanifold chaos increases the more the 
cyhndrical symmetry is broken and becomes strongest in 
perpendicular crossed fields. The interesting result of this 
Paper is that this assumption is not true and intramani- 
fold chaos becomes strongest, in our examples, at mutual 
field orientations around /3 « 60°. Effects of intramani- 
fold chaos are surprisingly weak at /3 = 90° as can be seen 
in Figs. 8e and lOe. Probably this cannot be explained 
by the existence of a discrete symmetry, i.e., the z-parity, 
TTz, in perpendicular crossed fields and further investiga- 
tions are necessary to understand the mechanism of the 
increase and decrease of intramanifold chaos. 

Implications are that some effects which have been dis- 
cussed recently for the crossed field atom such as Er- 
icson fluctuations of the continuous photo cross section 
above threshold |^ and Arnol'd diffusion of the three- 
dimensional classical motion might be much stronger 
at mutual field orientations around /3 ~ 60° than in per- 
pendicular crossed fields. It might be useful to revisit 
these problems and extend the previous investigations to 
arbitrary mutual orientations of the external fields. 

Up to now experimental investigation on atoms in com- 
bined electric and magnetic fields have concentrated on 
parallel and perpendicular field orientation. If indeed 
some interesting physical effects are much stronger for 
special mutual field orientations than in parallel or per- 
pendicular fields it is a challenge to observe these effects 
experimentally. 
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